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We present a supersymmetric analysis of the wave problem with a Demkov-Ostrovsky 
spherically symmetric class of focusing potentials at zero energy. Following a suggestion 
of Levai, we work in the so-called Ho = sector in order to obtain the superpartner 
(fermionic) potentials within Witten's supersymmetric procedure. General solutions of 
the superpotential for the known physical cases are given explicitly. 

We consider the wave/quantum problem at fixed null energy, (or zero-binding 
energy), [— |^V^ + UK{r)]^jj{r) = 0, with Demkov-Ostrovsky (DO) focusing poten- 
tials 

U (r) = ^ (1) 

^ (r/i?)2[(r/ii')---t-(r/i?)«]2' ^ ' 

where w>0, R>0, k>0 are constant parameters, and we introduced an energy 
scale £q — h? /2mR'^ for the potential part, following Daboul and NictoB. Of the C/^ 
class of potentials only the two cases to follow have been studied in the literature. 
For K = 1, one gets the wave problem for the Maxwell fish-eye (MF) lens. El whereas 
the K = 1/2 case has been used for the atomic Aufbau (AA) chart since it fulfills the 
Madelung rules of atomic energy ordering □. Demkov and Ostrovsky have shown 
that for the cases k = ki/k2, with fci and k2 integers, i) the classical trajectories of a 
zero-energy (i.e., zero velocity at infinity) particle close after k2 revolutions around 
the force centre, and ii) all the trajectories passing through a given point come to 
a focus after /c2/2 revolutions. 

The above wave equation can be written in the scaled variable p — r / R, (here- 
after the energy scale is to be understood), as follows 

' '^'^'^ " Uip) = 0. (2) 



5p2 pdp p2 p2(l-«)(l+p2K)2 

It is quite straightforward to solve the Sturm-Liouville problem, Eq. (2). Moreover, 
it can be turned into an eigenvalue problem of the DO coupling constant, w, and 
also written as a Laplace equation on the four-dimensional sphere EJ. The known 
results are the following. For w taking the quantized values, wn.k — {2k)'^[N + 
{2k)~^ — 1][N + (2k)^^], the regular, normalizable solutions read 

ll^Nlm{p) = RNl{p)Ylm{0, (f)) , (3a) 
r> _ ^(2^+l)/2K + l/2.,^ 
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where ^ 



N — n + {hi ^ — 1)Z, n = rir + ^ + 1, = 0, 1, 2, are the 'total', 



'principal' and 'radial' quantum numbers, I and m are the spherical harmonic num- 
bers, C^(^) are the Gegenbauer polynomials, i.e., the solutions of the corresponding 
ultraspherical equation (see Eqs. (7) and (8) below), and A/at; are the normalization 
constants. What one gets when the w parameter is made bigger and bigger is an 
increase of the degeneracy of the normalizable state at zero-energy, but only for the 
quantized values wm,h- The degree of degeneracy of such a group of states is N"^ , 
{N — 1,2,3...), similar to the electron energy levels in a Coulomb field. 

We pass now to the functions uni — pRni fulfilling the one-dimensional (half- 
line) radial equation 



H'u = 
with the effective potential 



u = 0, 



U, 



_ + _ {2kY[N + {2k)-^ -l\[N+ {2k)- 

eff - p2 - ^2(l-K)(l+p2K)2 



(4) 



(5) 



where we have already included supersymmetric superscripts. The functions uni 

'i(2i + l)/2K- 



are of the type /(p)C^itl^,''?^^^^^(^(p)), where f{p) reads 



f{p) = 



(^l ^ p2K)(2;+l)/2K ■ 



(6) 



and the Gegenbauer polynomials, C^, of degree p = rir and parameter q as given 
above, are the solutions of a second-order differential (ultraspherical) equation of 
the type 



dC_ 



Rp{OC{i) = 0, 



with 



Q{0 = 



21 + 2K+1 



and 



p{2q + p) 



1 



(7) 

(8.1) 
(8.2) 

(8.3) 



In Eq. (8.3), we emphasized the indexing of the R functions according to the var- 
ious sectors v = rir (0,1,2,...), which is a well-known characteristic of orthogonal 
polynomials □. 

We have now all the requisites for a Natanzon-type approach that we outline 
here following a recent discussion of Lcvai □. The method has been first used by 
Bhattacharjee and Sudarshan i and later by other authors, among whom Natanzon 
i is the best known due to his systematic treatment of hypergeometric cases. The 
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scheme deals with the fact that the solutions V' of any one-dimensional Schrodinger 
equation can be written as '4'{x) = f{x)F{z{x)), where f{x) is a function to be de- 
termined and directly related to the superpotential, while F{z) is a special function 
which satisfies a second-order differential equation of the form 



d^F ^, .dF 
dz^ dz 



1 ' Q{z)— + R{z)F{z)^Q. (9) 



In our case, Q{z) and R{z) corresponding to the Gegenbauer polynomials have been 
written above, while x = p and z{x) = £,{p). Then, the following equations can be 
readily obtained 

^ '%^Q{i{p)). (10) 



and 

/" U^ff 



mip)), (11) 



(ov icr 

where Ueff is given by Eq. (5). From Eq. (10) the function f{p) can be written as 
follows 

fip) - (e (P))-'/' exp [- j Qiapm] ■ (12) 

As suggested by Levai, one can define a ground state by the i?o(C) — sector, 
within which the Gegenbauer polynomials are Cq — 1 for any parameter q. In this 
simple case, from Eq. (11) one gets 

Ueff = +W\p) - — , (13) 

with W{p) — —-j^\nf{p). Eq. (13) is just the initial Riccati equation for the DO 
cases (and for any radial oscillator) in Witten's supersymmetric quantum mechanics 
0. Hereafter Eq. (13) will be called DORE. Since we actually know from Eq. (6) 
the function f{p) in the DO cases (one can check that Eq. (12) leads to the same 
function), a short calculation gives the DO superpotential 



/ 21 + 1 

P p{l + P^K) 

The effective DO superpartners in the Rq = Q sector can be written as follows 

dW^p) 
dp 

and 



W.ip) - z ~ TTT—^y (14) 



Ueff= ^ + W'Jp), (15a) 



''ff ' dp 
Thus, 



U-rf = +'^ + WUp). (156) 



_ ^(; + l) {2l + l){2l + 2K+l) 



4 H.C. Rosu 



and 

, (2Z + 1)(2/-2k-1) 2(2/ + 1) 

eff - p2(l-K)(l_^^2K)2 + p2(l+^2«) 

The factorizing operators read 



^ = + T^«, (17a) 
dp 

and 

^+ = -:^ + T^.. (17&) 

We have plotted V~j^ and V^^^ for some values of the parameters in Figs. 1 
and 2. From the plot of the DO fermionic potentials one can notice their repulsive 
nature. Consequently, the fermionic equation should be written in the continuum 

iZ+Mi = AA+ui = (--^ + i7+^)wi = fc^wi, k e (0, c3o). (18) 

It will be investigated elsewhere. Here we remark that in order to get the supersym- 
metric increment in the effective potential we used only the particular solution of 
the Riccati equation coming out from Eq. (11). On the other hand, it is well-known 
that the connection with the Gel'fand-Levitaji inverse scattering method requires 
the general solution of the Riccati equation eI. We construct it in the usual two 
steps as follows. Firstly, consider W = + VK(p) as another solution. Then by 
substituting in the DORE one gets 

dV 

— + 2VK(p)V = -1 . (19) 
dp 

This equation can be written as follows 



d_ 

dp 

with the solution 



Vexp ^2 j W{p)dp^ = -exp ^2 j W{p)dp 



(20) 



(21) 



V = - exp ^-2 j W{p)df^ ■ J exp ^2 J W{p)dp^ dp . 
Since we know that W{p) ^ ln/(p) we get J W{p)dp ^ ~ ln/(p). Thus 

V = -fip) I f-\p)dp . (22) 



The object of interest is now the integral of the inverse square of /, which reads 
explicitly 

^ dp . (23) 
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Let p'^ — tan(^)- Then the integral turns into the form 

2i + l 

^ / (24) 

K J (sina)(2'+'-+i)/« ' ^ ' 

and for the cases of physical interest, MF and AA, the formulas 2.515.1 and 2.515.2, 
respectively, in Gradshteyn and Ryzhik □ should be used to express it as a series. 

Thus, in the MF case, {n —\)., the integral Eq. (24) can be worked out into the 
series 

- -2m |(--)""^ + E [(2/-l)(2^-3)...(2^ + l-L)] ^--^"^ 

(25) 

while for the A A case, (k =1/2), the integral Eq. (24) reads 

24'+3 „^ ,,_.r ^ ^ [(4/ + l)(4;-l)...(4;-2m + 3)] i 4[(4Z + 1)!!] 



- cosa(csc2 a)(^'+i) h + V + ^Jl^/ - i)...(4/ - l^m + 3)j 1 4^4^ + iJ!!J « 

■^S- 2^ + 1^°'"^''" l'+;f^^(2csc2a)'"[(2/)(2/-l)...(2Z-m + l)]J+4-'[(2Z + l)!]'''*''''4^' 

(26) 

where a = 2 arctan p in the first formula and a = 2 arctan in the latter one. 

In the MF case, the radial factor /^(p) can be written trigonometrically as 
2~^'(sina)^' sin^(^) implying 

V _2cosa « / (2sin^a)"[^(^-l)...(^-m)] | 

- 27TT ^""''^ 2 ^ \ hi [(2? - mi - + 1 - 2m)] / ' 

In the AA case, the square radial factor is 2^"" (sin a)^^ sin''(^) and one can work 
out easily a formula for Vi 

2cosa^ 2.«,r, , ^ (^)'"[(4; + l)(4^-l)...(4;-2m + 3)] i , + 1)11]^-^^ , , 

2m'"" ^2Hl+E (20(2/-l)...(2;-m+l) J +^^1)1^1)]^ 2 ^ ' 

(28) 

We have in this way all the ingredients for the second step, which is to write down 
the general solution of DORE, containing a constant A-parameter in Eq. (22), see 
Refs. The general solution means Eq. (22) modified as follows 



Va,. = -mp) [^ + J f.iP )dp j . (29) 

In this case, the integral equation (24) will have the lower limit a and the upper 
one TT. Thus, Va.k can be calculated from our formulas as Va,k — Vk. — Xf^{p). 
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Figure Captions 

Fig. 1. Effective DO superpartners in the Rq = sector, U^jf{p), (a) and (b), 
respectively, in units of £o, for I = 2, and k = 1/2, 1, and 3/2. 

Fig. 2. Effective superpartners of the i?o = sector in £q units: (a), U~^j for 
1= 1, 5, 10, and (b), [/^^ for Z=6, 7, 8, in the case of Maxwell fish-eye, k = 1. We 
have plotted U'^ff in the region of the critical (inflexion) angular number, Zcr, that 
wc have found numerically to be Zcr=6.876 for /9cr=l-599. The critical / is the entry 
point toward a pocket (trapping) region of U^^^ for I > l^r- 

Note added on January 21st, 1996 

The mathematics in the last two pages of the paper (pp. 38, 39), though not 
wrong, might be considered as misleading with respect to the literature, since we 
obtained the general solution of the initial 'bosonic' Riccati equation. In this way, 
one can introduce the one-parameter family of fermionic potentials with the same 
bosonic superpartner. This is what we have done. However, for the usual con- 
nection with the Gcrfand-Lcvitan method, one should obtain the general solution 
of the 'fermionic' Riccati equation, and thus, the one-parameter family of bosonic 
potentials with the same fermionic superpartner. The task can be readily done on 
the base of the results at pp. 38, 39. 
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